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An algebraic method has been developed which allows one to engineer several energy levels in- 
cluding the low-energy subspace of interacting spin systems. By introducing ancillary qubits, this 
approach allows k-body interactions to be captured exactly using 2-body Hamiltonians. Our method 
works when all terms in the Hamiltonian share the same basis and has no dependence on pertur- 
bation theory or the associated large spectral gap. Our methods allow problem instance solutions 
to be embedded into the ground energy state of Ising spin systems. Adiabatic evolution might then 
be used to place a computational system into it's ground state. 



This work considers an important problem. Given 
a Hamiltonian comprised solely of 1-body and 2-body 
terms, from this Hamiltonian, and with the aid of an- 
cillary qubits, is it possible to construct the ground 
states of a Hamiltonian containing k-body terms with 
respect to a suitable subspace? In both the classical and 
quantum cases, this problem is particularly important 
when considering the physical complexity of interacting 
spin systems evolving into their lowest energy configura- 
tion [1, 2, 3, 4] or the equivalent computational task of 
determining the ground state [5, 6]. 

The ground state energy problem has long been con- 
sidered in the realm of classical complexity theory with 
well known results appearing in work such as [1, 5]. The 
extension to quantum complexity classes was prompted 
when Kitaev [6], inspired by ideas from Feynman [7], 
showed that the ground state energy problem of the 5- 
local (that is, 5-body) random field quantum spin model 
was complete for the quantum analogue of the class NP. 
Thus it was shown that 5-local Hamiltonian was 
QMA-complete and the quest to determine the complex- 
ity of various spin models began [10, 11, 12, 13, 14, 15, 
16, 23]. 

Ideas from the theory of quantum computation have 
also led to the use of ground state properties of quan- 
tum systems for computation [3, 8, 9]. This is known as 
the adiabatic model of quantum computation [3, 8] — in 
which a driving Hamiltonian is slowly replaced, most of- 
ten with a commuting Hamiltonian with a ground state 
spin configuration representing a problem instance solu- 
tion. 

At the heart of the construction of the QMA- 
completeness proofs lies the development of methods to 
engineer low-energy effective Hamiltonians, which ap- 
proximate k-body interactions, using at most 2-body 
terms [12, 13, 14]. To date, all known methods require 
the introduction of a large spectral gap, where the magni- 
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tude of the gap improves only an approximate low-energy 
effective Hamiltonian. It would be desirable if one could 
i.) remove the spectral gap dependence by capturing 
the low-energy effective subspace exactly and ii.) de- 
velop a systematic method to engineer multiple energy 
subspaces, including any ground state. 

The present paper addresses both of these problems. 
Somewhat surprisingly, it is possible to remove depen- 
dence on the large spectral gap by allowing the state of 
the ancillary mediator qubits (facilitating the coupling) 
to follow the state of the qubits being coupled. In appli- 
cation, care is taken to ensure that the active role of the 
mediator qubits is appropriate for any given application. 
In many cases, this new approach allows ground states of 
k-body interactions to be captured exactly using 2-body 
interactions; under the restriction that all terms in the 
Hamiltonian share the same basis. 

Structure: The remainder of this paper begins with a 
short introduction, followed by § II. which explains how 
the ground states of 3-body Hamiltonians can be used 
to embed any Boolean function (and for that matter, 
any switching circuit). § III reduces the 3- local Hamil- 
tonians used in § II to the case of 2-local Hamiltonians: 
In addition, we prove Theorem III.l, which states the 
existence of an efficient method to construct Hamiltoni- 
ans that simulate Boolean functions containing k- variable 
couphngs (i.e. Aa;2 A • • ■ AXk). In § IV we construct 2- 
body Hamiltonians that exactly capture the ground space 
of k-body Hamiltonians of the form Jai ® a2 ® ■ ■ • ® Ck- 
§ IV also contains a proof of Theorem IV. 1, which states 
the existence of a method to construct several energy sub- 
spaces of a given Hamiltonian — a necessity for certain 
applications. 

In addition to the main body of the present paper. 
Appendix A presents a proof of a tailored variant of 
the projection Lemma [9, 10, 12]. This is followed by 
Appendix B which explains Karnaugh maps — key to 
an algebraic reduction method relied on during several 
derivations. We make use of standard quantum comput- 
ing notation and background information [6, 10] as well 
as that for discrete functions and circuits [17, 18]. 



I. INTRODUCTION 

Let us represent an Ising spin with index i by the vari- 
able Si € {+1, — !}■ One could also represent variable Si 
in terms of binary variable Xi € {0, 1} as = 1 — 2xi, 
which we will denote as A single spin system can be 
acted on by linear combinations of operators taken from 
the set {Ijzbcr}, where the identity operator (1) can be 
scaled to ensure positive-semidefiniteness and the opera- 
tor a has eigenvectors |0) and |1) with respective eigen- 
values -1-1 and —1. The energy levels of the Hamiltonian 
operator ^(1 + en) (^(1 — (Ji) respectively) correspond- 
ing to the states |0) and |1) are 1 and (0 and 1 respec- 
tively). Addition of the operator i(l -I- ai) (^(1 — en)) 
adds an energy penalty to the state |0) (|1)) and can be 
thought of as negation (assignment) of variable Xi. 

In the case of two Ising spins, a complete basis of con- 
figurations are |00), |01), |10) and |11). Let us add scaled 
sums of a coupling term to our Hamiltonian: iataj. One 
can think of adding the operator ^(l — f^CTj) as a log- 
ical equality operation (i.e. the characteristic function 
true) and the operator 5(1 + cno'j) as a log- 
ical inequality operation (i.e. Xi <^#> Xj is true) between 
spins. For example, assume we act on a dual spin system 
with the Hamiltonian for inequality: the ground space 
is in span{|01), |10)}, so any vector that corresponds to 

two spin variables being equal (e.g. span{|ll), |00)} 
span{|x)|y)|a; = y,\/x,y £ {0,1}}) receives an energy 
penalty. 

We have shown how to set single spin variables, and 
how to apply equality and inequality operations between 
two spins. These operations, however, do not form a 
convenient logical system [24]. This will be done next, 
in § II and HI, by defining Hamiltonians with ground 
state spin configurations representing logical operations 
such as the AND (A) gate, the OR (V) gate, etc. We 
know that these dual arity operations require at least 
three spins as XiOxj = Z-i,. What we need is to find a 
way to set the low-energy subspace of three spins Si , sj 
and to be, for instance, the logical AND of the spins 
Si A Sj = Zi,. This assignment turns out to be possible 
working in the energy basis of a Hamiltonian equipped 
with a commuting local field and coupling term, such as 
an Ising Hamiltonian [25]: 

-falsing = ^ hi(7i + ^ JyCTjO-j. (1) 

Impressive demonstrations using qubits based on Joseph- 
son junctions [19, 20, 21] make an adiabatic [3, 8] real- 
ization of ground state logic gates using variants of the 
Hamiltonian (1) a foreseeable possibility. 

II. GROUND STATE SPIN LOGIC 

Consider some Hamiltonian H acting on a Hilbert 
space TL that is a sum of the vectors spanned by the 
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FIG. 1: Illustrating the mapping between circuits (with 
boolean variables Xi) and spins [si) for the example given 
in (3). One can use any number of methods to embed logical 
networks [17] into the ground space of Hamiltonians. 

subspace C and the orthogonal component of C written 
as thus Ti = C^ + C. The lowest eigenvalue of H will 

be denoted as X{H). Now let H^ ~ C) be defined 

as a projector onto £. Then H^iJH^ is the restriction 
of H to the subspace C — let us write this restriction as 
H\c- 

To develop the logic, consider the Hamiltonian TJpi-op 
such that Hproplc = and Hp^oplc^ > S {> 2||7Jin||) 
where Hin is a perturbation later used to set the circuits 
inputs, the norm || • || is the magnitude of the Hamiltoni- 
ans largest eivenvalue and S is the spectral gap between 
the and £ subspaces. We are faced with the task 
of ensuring that TJpropIc is a zero eigenspace when C 
spans the truth table of the logical operation of interest 
(e.g. C = span{|a;i)|a;2)|a;ina;2)|Va;i,X2 G {0,1}}). Let 
C be the low-energy subspace representing the truth ta- 
ble in the binary observablcs. Explicitly, in the case of 
logical AND, £. = span{|000), |010), |100), [111)} (ordered 
\xiX2)\Zi,), where z^, = X1AX2), which is a zero eigenspace 
of iJprop and = span{|001), [Oil), |101), |110)} (z^ ^ 
xi A X2) will be all eigenspaces of at least i5 [26]. 

One can add a perturbation, Hm, to set the circuits 
inputs. We will write this as a projector onto the n long 
binary bit vector x. This 1-local projector has the form: 

H.^=^^ix^)(x^i=Q)x:(i+(-i)""-o- 

Now upper bound ||-ffin|| (for all two input and single 
output gates [27]) as \\H\\ < 2. This implies that the 
spectral gap S is greater than 2. By noticing that k 

H\s,) = X\s,), H\si) = \\si) 
and {s,\H\si) + {si\H\s,) = 0, 

where \sj) G C and \sk) G C^, one recovers the strict 
equality, A(i?in|£) = A(iJ) (see Lemma II. 2). 

Using combinations of these ground state logic gates, 
we will perform computations. For example, write the 
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Hamiltonian with a low-energy subspace in 

span{|a;i.T2)|y*) | = a:i A X2, Vxi, X2 e {0, 1}} (2) 

as H/-,{xi,X2,yi,) and, with defined in (2), write the 
Hamiltonian with a low-energy subspace in 

span{|a;3?/^)|z^) | Aa;2,Va;3 £ {0,1}} 

as i?/\(a;3, j/^, z*). Then the proposition xi Aa;2 Vxs = 
is constructed as a sum of terms: 

/ . 

H = Hy{xi,x2,y^) + H^{x3,yi,,z^)+Hin (3) 

and the circuits input, iJjn, is yet to be defined. The 
qubit with label z^, is now equal to cci A 0:2 V X3 and is 
a temporary variable that is equal to xi A X2, as seen in 
Table II. 

A small perturbation, Hin, can be added to set any 
of the qubits to specified values. For example, to set 
the input as xi = 1, 2:2 = and 2:3 = one adds the 
perturbation i7i„ |0)(0|i + |1)(1|2 + |1)(1|3. If, alter- 
natively, we were to let iJjn = |0)(0|*, which acts on the 
circuits output z*, then the low-energy subspace would 
be spanned by all vectors where the output z^, is |1) [28]. 
As seen from Table II, this subspace is in 

span{|001)|l}|0), |011)|1)|0), |101)|1)|0), 
|110)|1)|1),|111)|1)|1)}, 

where we adhere to the ordering \xiX2X3)\zi,)\yi,) . If in- 
stead we were to add the perturbation Hi^ to the qubit 
labeled jy^,), the ground space would be spanned by 
{|110)|1)|1),|111)|1)|1)}. ^ 

To complete our reduction, the 3-local Hamiltonians, 
just described, will be reduced in the next section to 2- 
local Hamiltonians. Before continuing to our 2-local re- 
duction, let us state Lemma II. 2 and Theorem II. 1 — 
the proof of which is implied by the results of this sec- 
tion. Here we choose a finite set f2 of one-output Boolean 
functions as basis. Then, an fJ-circuit works for a fixed 
number of Boolean input variables and consists of a finite 
number of gates, where each gate is defined by it's type 
taken from f2. (For additional background information 
on boolean functions and switching circuits see the freely 
available standard reference [17].) 

Theorem II. 1. Let f be a switching function given as 
the map f : {0, 1}'^' {0, 1}™ for finite k andm. Now let 
there be an asynchronous Q- circuit computing f. Then 
there exists an fl- circuit embedding into the ground space 
of a 3-local Hamiltonian, H3, such that: i.) The norm 
of the Hamiltonian ||i?3|| is constant and, in particular 
independent of the size of f , the ^.-circuit, as well as k 
and m. ii.) The ^.-circuit embedding is upper bounded 
by a number of qubits 0{1)- reducible to the number of 
classical gates required on the same lattice. 

An important technical tool used in our construction is 



a variant of the projection Lemma [9, 10, 12] — proven in 
Appendix A. Let us denote H as a Hilbert space of inter- 
est and let Hi be some Hamiltonian. Consider a subspace 
C € Ti. such that a Hamiltonian H2 has the property that 
is a eigenspace and C-^ is an eigenspaee of at least 
6 (> 2||i7i|[). Consider the Hamiltonian H = Hi + H2, 
the projection lemma says that the lowest eigenvalue of 
H, X{H), is the lowest eigenvalue of Hi restricted to the 
subspace C — that is X{Hi\c)- Thus, by adding H2 one 
adds a penalty (proportional to S) to any vector in 
To state the Projection Lemma (Strict Equality) we: 

Lemma II. 2. Let H = H1+H2 be the sum of two Hamil- 
tonians operating on some Hilbert space TL = C -\- . 
Denote C = span{|sj)|Vj} and = span{|s^)Vfc} 
for finite j,k. Consider the restriction H2\c = and 
H2\c^ > S{> 2\\Hi\\). Then, if\/j,k H\sj) = X\sj) (V/c 
^14) = Msi)), {sj\H\si) + {si\H\s,) = the following 
equality holds: X{H) = X{Hi\c). 



III. THE 2-LOCAL REDUCTION 

The main result of this section can be found in Ta- 
ble III. To develop this table we used the algebra of 
multi-linear forms [18] and the Karnaugh map method 
from discrete mathematics [22] — which we review in 
Appendix B. 

We consider multi-linear forms that are maps / from 
the Booleans numbers to the reals, where the inputs and 
outputs are of finite size. For instance, the multi-linear 
form for AND (OR) is simply = a;iAa;2 (/v = Xi-\-X2 — 
2xi A X2). Hence, one can express the Boolean equation 
f — xi A X2\/ X3 with the polynomial f = xi A X2 -\- x^ — 
xi A a:;2 A 2:3. Let us first write the vector of integers: 

C"^ = (co,Ci,C2,C3,C4,C5,C6,C7,), (4) 

representing the outputs of a multi-linear function / 
over the three Boolean input arguments xi , a:;2 and X3 . 
We wish to construct a canonical representation for any 
multi-linear function of three variables in terms of the 
vector c from (4). We will represent the negation of the 
variable (or using the notational equivalent -ix) 

and canonically expand (4) as a sum of products: 

f{xi,X2,X3) = Co ■ X1X2X3 -\- Ci ■ X1X2X3 -\- C2 ■ X1X2X3 

+C3 ■ X1X2X3 -\- C4 • X1X2X3 -\- C5 ■ X1X2X3 
-\-ce ■ X1X2X3 -\- C7 ■ X1X2X3. (5) 

This expansion (5) forms a basis for the space of 3- 
variable Hamiltonians, but to realize any of the eight 
terms requires 3-body couplings. This motivates us to 
write a second canonical expansion, found from a change 
of variables in (5) and by expanding each term into it's 
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positive polarity form: 



f{xi,X2,X3) = aa + ai ■ xi + a2 ■ X2 + as ■ X3 
+ 04 • X1X2 + as • X1X3 + ae ■ X2X3 
+ 07 • X1X2X3. (6) 

This equation (6) also forms a basis for the space of re- 
alizable Hamiltonians of 3-spins. In this suggestive form, 
however, we can truncate (6) past 2'"' order and con- 
sider the subclass of Hamiltonians that can be realized 
by setting 07 = 0. 

Out of the 16 possible functions of 2-input and 1- 
output variable, it can be proven that only two are not 
realizable using 3-spins. These are the 2-local penalty 
Hamihonians for XOR (e) and EQV (0) [29], which are 
each possible to realize by adding a single mediator qubit 
(as seen in Tabic III). 

We will explain in detail how the positive-semidefinite 
AND penalty Hamiltonian, H/^, is derived. We anticipate 
that the details of our approach will aid others faced with 
Hamiltonian constructions. Let C be the null space of 
and let all higher eigenspaces be given as C-^. The 
penalty Hamiltonian has a null space, C, spanned by the 
vectors {|xiX2}|z*}|zi = xi A 2:2, Vxi, a;2 € {0,1}}. De- 
note S as an energy penalty applied to any vector com- 
ponent in . Our goal is to develop a Hamiltonian that 
adds a penalty of at least S to any vector that does not 
satisfy the truth table of the AND gate — that is, we want 
to add an energy penalty to any vector with a component 
that lies in C^. 

In order to make the penalty quadratic, one first con- 
structs the Karnaugh map illustrated in Fig. 2 c.) for the 
case xi Ax2 = z*- This is done by examining Table I. In 
the right most column, all possible assignments for the 
variables xi, X2 and Zi, are shown. The Karnaugh map 
is constructed by examining the second column. When- 
ever the variable is not equal to the AND of the vari- 
ables xi and X2^ a penalty of at least 5 must be applied, 
which ensures that vectors in the ground space satisfy 
1x2)1x1 A X2). Any vector that must receive an en- 
ergy penalty of 5 is depicted in the Karnaugh map with 
a dot (■). 

Begin by noticing that any vector associated with cube 
number 4 must receive an energy penalty, so the 1-local 
field corresponding to the qubit with label z^, must be at 
least 8 — adding the term piz^, to the Hamiltonian, with 
the constraint pi > S. Cube 3 must also receive an energy 
penalty of at least 6, adding the term P2X1 A X2 to the 
Hamiltonian H/^. With both penalties applied, vectors 
corresponding to cube 7 must be brought back to the 
null space — accomplished by subtracting the quadratic 
energy rewards riz^Axi and r2Zi,t\X2 from iJ^^. A system 
of equations for the Hamiltonian H^(xi, X2, z^,) — 



piz^ + P2X1 hX2- r2Z* A Xi - r2Z^, A Xi 



(7) 



can be solved to set the rewards (r's) and the penalties 
(p's). This system is derived from the fact that the term 
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FIG. 2: Karnaugh maps: a.) 2-local (positive polarity) inter- 
actions circled (e.g. qiXiX2+q2XiAZi,+q3X2Af*). b.) Linear 
(positive polarity) fields circled (e.g. lixi + I2X2 The 
interactions in cubes a.) and b.) form a basis for the space of 
realizable (< 3 qubit) positive-definite logical gadget Hamil- 
tonians expressible as: H{xi,X2, z^^) — ko + hxi + I2X2 + 
hx3 + qixi /\X2 + q2Xi /\Zt, + q-iX2 /\Zi,, where Vi, fco, /i, qi > 0. 
c.) A Karnaugh map illustrating (with ovals) the linear and 
quadratic terms needed to set the null space of the Hamil- 
tonian (8) to be in span{|a;i2;2)|j/*)|i/+ = xi A X2,^x-i,X2 £ 
{0,1}}. 



X1X2X3, corresponding to cube 7, must have zero energy: 
= pi 4-p2 — ''I — ''2 and is subject to the conditions that 
P1jP2 > S and |ri + r2\ > Pi- For convenience, let 6 = 1 
and then determine values for the coefficients in (7) and 
thus derive the 2-body Hamiltonian (for AND): 

H/\{xi,X2, 2*) = 3z* + xi Ax2~ 2z* Axi — 2z* A X2- (8) 

If one desires to invert an input variable, she simply 
applies the transform: Xi —^ (1 — Xi). For example, the 
Hamiltonian applying the penalty Hf^{xi, X2, z*) is: 

3z^ + (1 - xi) A 2:2 - 2z* A (1 - xi) - 2z* A X2- (9) 

To write this Hamiltonian in terms of spin variables, 
first change each variable, Xi, to its (matrix) operator 
form by the replacement Xi — > |0)(0|i. The change to 
spin variables is then accomplished by the replacement: 
|0)(0|i ^(1 — (Ti). After these substitutions one arrives 
at the Hamiltonian (10) which is isomorphic (~) to (9). 



H^{xi,X2,Zi,) ~ iJA(-Sl,S2,S*) 
1 



(10) 



-(3 + CTl - 0-2 + 2(7* - (7l(T2 + 2(71 tr* - 2(72(7*) 

4 

We now have the necessary machinery in place to state 
two theorems (III.l and HI. 2). In the first, we are con- 
cerned with a situation that arises in several applications. 
That is, one often needs to couple three Boolean variables 
(AND product), as Xi A 0:2 A 0:3, using only 2-local Hamil- 
tonians. From our reduction, it is possible to efficiently 
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TABLE I: Left column: possible assignments of the variables 
Xi, X2 and 2*. Center column: illustrates the variable as- 
signments that must receive an energy penalty > S. Right 
column: truth table for H/\{xi,X2, z*) = 3z* + xi AX2 — 2z* A 
xi — 2zi, l\X2^ which has a null space £ e span{|a::ij:2)|2:*)|2:* — 
Xi A X2,'^Xi,X2 G {0, 1}}. 

construct any k-local product term, xi A X2 A ■ ■ ■ A Xk, 
of this type. We prove this in Theorem III.l. Wc then 
present Theorem III. 2, which is a 2-local variant of The- 
orem ILl — the proof of which follows directly from the 
results of this section. 

Theorem III.l. Let /j. be a k-local multi-linear from 
and let there be a Hamiltonian H^. acting on the Hilhert 
space Tik such that fk — Hk- Then there exists a 2- 
local multi-linear form, f2, and corresponding Hamil- 
tonian, H2, acting on the Hilbert space Ti.2 (where 
Ti-k C 0.2), with the same low-energy subspace of H2 in 
span{|x)|y}|y = fk{x),Wx e {0,1}", Vy G {0,1}™} C H. 
The number of mediator qubits required to realize H2 is 
upper bounded by 0{size{fk)). In addition, the spectral 
gap of H2 is bounded by the spectral gap of H^. . 

Proof. To construct such a Hamiltonian, we will employ 
an inductive argument and consider a single (out of w) 
k-local term, hk = xi A X2 A ■ ■ ■ A Xk, that couples k > 3 
Boolean variables. We will now show the existence of 
a 2-local reduction requiring (k — 2) mediator qubits to 
embed hk into the ground state of a 2-local Hamiltonian. 
Consider the 2-local coupling A X3 and add the Hamil- 
tonian that forces an energy penalty whenever z^, is not 
the Boolean AND of the variables xi and X2- The 2-local 
Hamiltonian is written as 

H/s,{xi,X2,Z^,) + z^, AX3 ~ i(4 - 0-1-0-2 + 

3ct* + 0-3 + 0-10-2 - 2ctio-* - 20-20"* + cr^o-a), 

where is found in Table IH, and z^, is a tempo- 
rary variable. In words, the variable is coupled to 
xg and the penalty, H^, forces to be equal to the 
Boolean product of xi and X2 — thereby creating the 
desired coupling with respect to the subspace spanned 
by \x1X2X3), Mi G {1, 2, 3}, Xi G {0, 1}. For a k-local term 
xi A X2 A • • • A Xfc, this procedure is recursively repeated 
k — 2 times. The reduction requires w{k — 2) qubits to 
capture the low- lying eigenspace of Hk with H2. □ 



Theorem III. 2. Let f be a switching function with a 
fixed number of inputs k and outputs m. Let there be 
an asynchronous ^l- circuit computing f over the basis 
{^,®,A}. There exists an il-circuit embedding into the 
ground state of a 2-local Hamiltonian, H2, such that: i.) 
The norm of the Hamiltonian ||i?2|| is constant and, in 
particular independent of the size of f , the ^-circuit, k as 
well as m. ii.) The ^.-circuit embedding is upper bounded 
by a number of qubits 0{k) -reducible to the number of 
classical gates required on the same lattice. 



IV. A NOVEL 3-LOCAL GADGET 

Wc are concerned with constructing the ground state 
of the operator Joi (2) 0-2 (8) 03 — which is a different task 
than coupling (that is, the AND product) three Boolean 
variables Xi A X2 A 0:3. Without loss of generality, let us 
consider construction of the target Hamiltonian 

^^targct + Jai (g) 0-2 <8) 0-3, (11) 

where Y is diagonal in the a basis. Wc will write the 
spectrum of o^i ® 0-2 C3i 03 , in canonical (Boolean counting) 
order, as {1,-1,-1,1,-1,1,1,-1} [30]. Now the low- 
energy, A(o-i ® 0-2 80-3) = —1, eigenspace is given as 

C = span{|001), |010), |100), |111)} 

and the high-energy, -1-1, eigenspace as 

^ span{|000), |011), |101), |110)}. 

Over the complex field, the tensor product (o) of two 
elements is simply their complex multiplicative (■) prod- 
uct. With respect to the canonical order, the spin 
variables for this operator (11) form the product z^ = 
si • S2 • S3, where Vz, Si G {+1, —1}, and so we consider 
the group Homomorphism { — 1,-|-1, ■} 1— s- {1,0,©}, where 
© denotes modulo 2 sum (XOR); whence 

z* = a;i © 2:2 © X3, Mxi,X2,X3 G {0, 1}^. 

In what follows, we will present a general framework to 
construct the ground state of any operator in the a basis 
and apply this approach to produce a 3-local gadget re- 
quiring three mediator qubits. We will then focus our at- 
tention on optimization of this new 3-local gadget, which 
is shown to be possible to realize using only two mediator 
qubits. 

Let us state an overview of our approach. To capture 
both the low- and high-energy spectrum, while preserv- 
ing the spectral gap, one will first write down a penalty 
Hamiltonian for the 3-variable function z*, which acts 
on the Hilbert space TC. This function, z^, outputs log- 
ical for any input vector in C, and for all vectors in 
the function outputs logical 1. We will next add a 
small perturbation to the output z^, — thereby breaking 
the low-energy degeneracy and allowing us to capture the 
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spectrum of (11) exactly, with respect to the subspacc 
C + £^ = span{|a;iX2X3)|Va;i, X2, X3 G {0, 1}} C H. 

3-local gadgets with 3 mediator qubits: From 
Table III we know that each XOR function requires 
an extra qubit, and so three mediator qubits are re- 
quired to create the desired coupling. Let us write the 
Hamiltonian that applies the XOR penalty to the vari- 
ables x\ and X2 as HQ{xi,X2,yn,mi) and the Hamilto- 
nian that applies the XOR penalty to the variables X3 
and Tji, as _ff (0:3, y^, z^, 7/12). Now order the variables 
as \xiX2X3)\Zi,)\yi,mim2), where mi and TO2 are medi- 
ator qubits and is a temporary variable that is not 
read. To split the spectrum into it's respective low- 
energy {£) and high-energy (£^) subspaces we add the 
perturbation V = J(|0)(0| — |1)(1|), which acts on the 
qubit z^. This allows one to construct the Hamiltonian 
(11), with the desired spectrum since the commutator 
[Y, Jai (g) (72 8) 0-3] = shows that Y only adds energy 
shifts and not level mixing — see Lemma II. 2. 

3-local gadgets v^rith 2 mediator qubits: Let us 
present an alternative approach to realizing a 3-local 
gadget which requires only two mediator qubits. To 
construct the gadget Hamiltonian, consider the 2-local 
coupling Z4S3 and add the Hamiltonian that forces a 
penalty whenever 24 is not equal to the EQV of vari- 
ables X2 and 3:3. The 2-local Hamiltonian is written as 
Hq{xi,X2, z'^) + -2453, where TJ® is found in Table HI, 
and z^ is a temporary variable. The Hamiltonian (12) 
captures the desired spectrum for S > 2\J\. 

H = -(4 + (72(73 + ((72 +(73)^4 + (12) 

2(1 - (72 - (73 - CT4)ct5 — (72 — (73 - CT4) + JoiOi . 

S1Z4 

The ground space of the Hamiltonian (12) is given as 

C = span{|001)|00), |010)|11), |100)|11), |111)|01)} 

and the first excited space as 

= span{|000)|01), |100)|00), |110)|00), |110)|10)}, 

where the qubits are in ascending order: qubit 4 repre- 
sents the Boolean EQV of qubits 2 and 3, while qubit 5 
is the mediator qubit needed to construct the EQV func- 
tion. 

We will now state then prove Theorem IV. 1 which al- 
lows one to construct, not only the ground state, but 
several levels of the low-lying energy subspace of k-body 
interactions using only 2-body Hamiltonians, formally we 

Theorem IV. 1. Let Hk be a k-local Hamiltonian diag- 
onal in any basis a and let this Hamiltonian act on the 
Hilbert space Tik- Assert that Hk has bounded norm, and 
let the strictly increasing list E2, ■ ■ ■ , denote the 



eigenenergies of formed by combing degeneracies, and 
label the corresponding eigenspaces as {Ci,C2, - • ■ T^k}, 
respectively. Then there exists a 2-local Hamiltonian, 
H2, with a low-lying spectrum isomorphic to that of Hk- 
Moreover, H2 is equivalent to Hk with respect to a sub- 
space spanned by {Ci, C2, ■ ■ ■ ,^k}- In particular, there 
exists a 2-local reduction capturing the k energy subspaces 
{Ci, C2, ■ • • , Ck} in the low-energy subspace Hk- 

Proof. Let us review the general method to construct 
ground states. First, determine C, the low-energy sub- 
space, and let Eg denote the ground state energy. One 
will next write a function, = /(xi, a;2, . . . , a;„), that 
outputs for all input vectors in C, and for all other vec- 
tors the function will output 1. The ground state will 
be realized with respect to a subspace spanned by the 
qubits labeled |xi)|x2), . . . , |a:,i)- To capture the desired 
ground space, a perturbation {V = Eg\Q){0\) is added, 
which only acts on the qubit z-^. Assume that we are 
instead interested in capturing several energy subspaces, 
with energies {Ei,E2,--- ,Ek}, and let us label these 
spaces as {Ci,C2,--- j-Cfc}, respectively. We will con- 
struct a function with k outputs, and repeat the process 
outlined above — this time acting on each respective j*'' 
function with the perturbation V = ^ 



V. CONCLUSION 

In this study, we have adapted a range of classical 
algebraic reduction methods to the construction of the 
low-lying energy subspace of k-local Hamiltonians us- 
ing 2-local Hamiltonians. Our methods do not rely on 
perturbation theory or the associated large spectral gap. 
We have examined explicit constructions of various use- 
ful k-local to 2-local conversion Hamiltonians — includ- 
ing both those needed to embed logical functions as well 
as couple spin variables. We have found constructions 
of these Hamiltonians which are optimal in the num- 
ber of introduced ancillary qubits. For ease of reference, 
our results are summarized in Table HI and Table IV. 
In Theorem IV. 1 we presented a novel method to con- 
struct several levels, including the ground state, of the 
low-lying energy subspace of k-body interactions using 
2-body Hamiltonians. Our methods have several appli- 
cations in adiabatic quantum algorithm design and quan- 
tum complexity theory. 
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APPENDIX A: PROJECTION LEMMA 

We will now prove Lemma IL2 which is discussed on 
page 3 in § IL Consider first the case that X{H) < 
^iHilc). Denote by [77) G £ the minimizing eigenvec- 
tor of Hi\c with eigenvalue X{Hi\c). Since H2\r]) = 0, 

(rylHjr,) = (r;|i/i|77) + {ti\H,\i^) - \{H,\c). 

Now consider actually minimizing over all vectors |C) of 
unit length: 

uAn {m\0}<m\v)^KHi\c). 
]C>e C+c^ 

proving a R.H.S. To show the lower bound on \{H) write 
any unit vector \v) S H = £ + £^ as \v) = a\s) + (3\s-^) 
where |s) {\s^)) is in C (C^), a, (3 e R, a, (3 > and 
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a^+l3^ = 1. So 
X{H) = \{Hi+H2)>a^{s\Hi\s)+a(3{{s\Hi\s^) + 

For real Hi, and 

a(3{{s\Hi\s^) + {s^\Hi\s)) = 2al3{s\Hi\s^). 

However, \s) and \s^) are eigenstates of i?i and (s|s^) = 
0, hence: 

X{Hi + H2) > \{Hi\c) + /3''{S - 2\\Hi\\) 
is minimized with /3 = so the projection lemma becomes 
HHilc) < X{H) < X{Hi\c) X{H) = X{Hi\c). □ 

APPENDIX B: KARNAUGH MAPS 

The Karnaugh map is a tool to facilitate the algebraic 
reduction of Boolean functions. We made use of this tool 
in § III during explanation of the specific details required 
to construct Tables III and IV. Many excellent texts and 
online tutorials cover the use of Karnaugh maps such 
as the wikipedia entry (http : //en . wikipedia . org) , the 
articles linked to therein as well as the straight forward 
reference [22]. This Appendix briefly introduces these 
maps to make the present paper self contained. 

Karnaugh maps (see Fig 2 for three examples), or more 
compactly K-maps, are organized so that the truth table 
of a given equation, such as a Boolean equation (/ : B" — > 
B) or multi-linear form (/ : B" M), is arranged in a 
grid form and between any two adjacent boxes only one 
domain variable can change value. 

This ordering results as the rows and columns are or- 
dered according to Gray code — a binary numeral system 
where two successive values differ in only one digit. For 
example, the 4-bit Gray code is given as: 

{0000, 0001, 0011, 0010, 0110, 0111, 0101, 0100, 1100, 
1101, 1111, 1110, 1010, 1011, 1001, 1000}. 

By arranging the truth table of a given function in this 
way, a K-map can be used to derive a minimized function. 



To use a K-map to minimize a Boolean function one 
covers the I's on the map by rectangular coverings con- 
taining a number of boxes equal to a power of 2. For ex- 
ample, one could circle a map of size 2" for any constant 
function / = 1. Fig 2 a.) and b.) contain three circles 
each — all of 2 and 4 boxes respectively. After the I's 
are covered, a term in a sum of products expression [17] 
is produced by finding the variables that do not change 
throughout the entire covering, and taking a 1 to mean 
that variable (xi) and a as its negation (xi). Doing this 
for every covering yields a function which matches the 
truth table. 

For instance consider Fig 2 a.) and b.). Here the boxes 
contain simply labels representing the decimal value of 
the corresponding Gray code ordering. The circling in 
Fig 2 a.) would correspond to the truth vector (ordered 
z^, xi then X2) 

(0,0,0,1,0,1,1,1)^. (Bl) 

The cubes 3 and 7 circled in Fig 2 correspond to the 
sum of products term XiX2- Likewise (5,7) corresponds 
to 2^X2 and finally (7,6) corresponds to Zi,xi. The sum 
of products representation of (Bl) is simply 

f{Zi,, Xi,X2) = X1X2 V Zi,X2 V 2^X1. 

Let us repeat the same procedure for Fig 2 b.) by again 
assuming the circled cubes correspond to I's in the func- 
tions truth table. In this case one finds 2* for the cir- 
cling of cubes ladled (4,5,7,6), X2 for (1,3,5,7) and xi for 
(3,2,7,6) resulting in the function 

/(2*, Xi,X2) = Xi V 2^ V X2. 

Our use of K-maps in § HI allows one to visualize cube 
groups (variable products) that are at most 2-local in size 
— the highest order terms realizable with 2-local Hamil- 
tonians. In addition, K-maps help reduce the number of 
simultaneous equations that, as seen in § III, must be 
solved — see (7) and (8). The Karnaugh maps shown in 
Fig. 2 a.) and b.) illustrate groupings for quadratic and 
linear interactions, respectively corresponding to 2-body 
terms and 1-local fields. In § HI, this observation al- 
lowed us to derive 2-local Hamiltonians and prove which 
Hamiltonians are not possible to construct given specific 
numbers of mediator qubits. 
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\X1X2X-J,) 




\y*) 


jOOO) 


|0) 


|0) 


1001) 


ll) 


|0) 


1010) 


|0) 


|0) 


1011) 


ll) 


|0) 


1100) 


|0) 


|0) 


1101) 


ll) 


|0) 


1110) 


ll) 


|i> 


1111) 


ll) 





TABLE II: Ground state truth table generated for the Hamiltonian (3). The function output, z*, is equal to a;i A 2:2 V xa. It is 
instructive to think of the variable y* as a coupler that follows the variables xi and X2 as = si A 0:2. 



function 


2-local Hamiltonian H(a;i, 2:2, 2*) = 


ground state (ordered: a;i) |a;2) |2*)) 


= 2* 

1 = 2* 


1(1+ as) 


span{|a;ia;2)|0)|Va;i,a;2 £ {0,1}} 
span{la:ia;2)|l)|Va;i, a;2 £ {0,1}} 


Xl f\X2 = 2* 
Xl /\X2 = 2* 
Xl /\X2 = 2* 
Xi /\X2 = 2* 


1(3 + (Tl + (T2 — 2(T* + CriCr2 — 2(71(7* — 2(72(7*) 
1(3 + (7l — (72 + 2(7* — (71(72 + 2(7l(7* — 2(72(7*) 
i(3 — CTl — (72 + 2(7* + (71(72 — 2(7l (7* — 2(72(7*) 
■|(3 — (7l + (72 + 2(7* — (71(72 — 2(7l (7* + 2(72(7*) 


span{l001), |010), |100), |110)} 
span{lOOO), |011), |100), |110)} 
span{!000), |010), llOO), llll)} 
span{lOOO), |010), |101), lllO)} 


XlW X2 = 2* 
Xl\/ X2 = 2* 
Xl V X2 = 2* 
XiW X2 = 2* 


^(4 + (7l + (72 — 2(7* + 2(71(72 — 3(71(7* — 3(72(7*) 
i(4 + (7l — (72 — 2(7* — 2(71(72 — 3(71(7* + 3(72(7*) 
^(4 — (7l — (72 + 2(7* + 2(71(72 — 3(71(7* — 3(72(7*) 
i(4 — (7l + (72 — 2(7* — 2(71(72 + 3(71(7* — 3(72(7*) 


span{lOOO), |011), |101), |111)} 
span{|001), |010), |101), |111)} 
span{i001), |011), |101), lllO)} 
span{|001), |011), |100), llll)} 


Xl 2* 
2:2 2* 
<ti> 2* 
X2 <S> 2* 


i(l+(7ia3) 
i(l -(72(73) 
i(l -(71(73) 
i(l +(72(73) 


span{|0a;2l), |la::20)|Vs:2 £ {0, 1}} 
span{|2;i00), |2:ill)|Vii G {0, 1}} 
span{|0a:20), lla-2l)|Va;2 G {0, 1}} 
span{|xi01), |xilO)|Vxi G {0, 1}} 


Xl ® a;2 = 2* 

XlQX2 = 2* 


4 + aiCT2 + ((7l + (72)(7* + 2(1 — CTl — CT2 — CT*)ct4 — CT2 — CT* — CT4 
4 — CT1CT2 + (cti — CT2)ct* + 2(1 — CTl + CT2 — CT*)ct4 + CT2 — CT* — CT4 


span{lOOOO), lOlU), |1011), lllOl)} 
span{|0100), lOOU), |1111), llOOl)} 



TABLE HI: Logical gadgets (§ HI): The span of the zero energy ground space (£) of these gadget Hamiltonians represent the 
truth table of a given switching function in the spin variables (as, for instance, the AND function: C = span{|a;i2;2) |2*) Iz* = 
Xl A 2:2, Vsi, X2 G {0, 1}}). This table includes all 16 = 2^ possible boolean functions with n = 2 inputs. 



3-local coupling 


2-local Hamiltonian 


Jxi A a-2 A 13 




CTl — CT2 + 3ct* + CT3 + CT1CT2 — 2ctict, — 2ct2CT* 


+ Jct*CT3) 


Jcti ® CT2 ® CT3 


|(4 + CT2CT3 + 


(ct2 + CT3)ct4 + 2(1 - CT2 - CT3 - CT4)ct5 - CT2 - 


CT3 — CT4) + JCT1CT4 



TABLE IV: 3-local gadgets. Top (§ HI): Hamiltonian with a low-energy subspace that couples three Boolean variables. The 
state of the mediator qubit ct* is a function (the AND) of qubits 1 and 2. Bottom (§ IV): Hamiltonian with low-energy subspace 
that couples three spin variables for 5 > 2|J|. The ground space, £ = span{|001)|00), |010)|11), 1100)|11), |111)|01)} and the 
first excited space, = span{|000)|01), |100)|00), jllO)|00), 1110)110)} — the qubits are in ascending order: qubit 4 represents 
the Boolean EQV of qubits 2 and 3, while qubit 5 is the mediator qubit needed to construct the EQV function. 



